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Topological doping and the stability of stripe phases. 
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We analyze the properties of a general Ginzburg-Landau free energy with competing order pa- 
rameters, long-range interactions, and global constraints {e.g., a fixed value of a total "charge") to 
address the physics of stripe phases in underdoped high-T c and related materials. For a local free 
energy limited to quadratic terms of the gradient expansion, only uniform or phase-separated config- 
urations are thermodynamically stable. "Stripe" or other non-uniform phases can be stabilized by 
long-range forces, but can only have non-topological (in-phase) domain walls where the components 
of the antiferromagnetic order parameter never change sign, and the periods of charge and spin 
density waves coincide. The antiphase domain walls observed experimentally require physics on an 
intermediate lengthscale, and they are absent from a model that involves only long-distance physics. 
Dense stripe phases can be stable even in the absence of long-range forces, but domain walls always 
attract at large distances, i.e., there is a ubiquitous tendency to phase separation at small doping. 
The implications for the phase diagram of underdoped cuprates are discussed. 
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INTRODUCTION 

One of the fundamental issues in the theory of highly 
correlated solids is the nature of the ground-state phases 
produced when a small concentration, x, of "doped holes" 
is introduced into a Mott insulator, Mfiicularly an an- 
tiferromagnct. It is now established^ □ that, at small 
enough x and in the absence of long-range Coulomb in- 
teractions, a doped antiferromagnet generally phase sep- 
arates into a hole- rich and a hole-free phase, i.e. the 
antiferromagnetic state is destroyed via a first order 
phase transition. In the presence of weak, long-range 
Coulomb interactions, that frustrate this local tendency 
to phase separation, the two-phase region is replaced by 
states pWjhich are inhomogeneous on intermediate length 
scalesfTB and especially "stripe phases," which have nqsa. 
been observed in a wide variety of oxide materials.Lrtll 
In various quasi two-dimensional cuprate high tempera- 
ture superconductorp-and the isostructural nickelates the 
stripes are observedli3 to be "topological," in the sense 
that the charge is concentrated along one-dimensional 
"rivers" which are at the same time antiphase domain 
walls in the antiferromagnetic order. In the nearly 
bic manganate colossal magnetoresistance materials 
the "stripes" are two dimensional sheets of charge which 
are non-topological. (In some sense, each ahjeet can be 
thought of as a dimer of topological stripes .Bo) 

Here we study the properties of a general Ginzburg- 
Landau free energy with competing order parameters, 
long-range interactions, and global constraints [e.g., a 
fixed value of a total "charge" , as defined in Eq. (g)] to 
address the physics of inhomogeneous ("stripe") phases. 
Specifically, a stripe phase is a unidirectional density 



wave which, in the case of a doped antiferromagnet, con- 
sists of a coupled spin-density wave (SDW) and charge- 
density wave (CDW). At very dilute doping, a stripe 
phase consists of an ordered array of far-separated self- 
localized structures, or individual stripes. At moderate 
doping levels, where the spacing between stripes is com- 
parable to their width, the structures are best described 
as nearly harmonic density waves. 

Zachar and two of ustil have considered the density 
wave limit of a Landau theory of coupled CDW and SDW 
order, each with a fixed wave vector q, near a transition 
to a disordered state, that occurs as the temperature or 
doping are varied. The existence of a cubic term in the 
Landau free energy coupling these two order parameters 
drives the period of the SDW to be twice that of the 
CDW, and the absence of any net AF ordering is equiv- 
alent to the statement that the stripes are topological. 
By contrast, as shown in Appendix |Aj, the same sort of 
term in the Landau theory of the transition between a ho- 
mogeneous ordered antiferromagnetic phase and a stripe 
ordered phase produces a state in which the antiferro- 
magnetic magnetization does not change its sign between 
the domains, i.e. the stripes are non-topological. 

To elucidate the circumstances in which arrays of 
stripes can be thermodynamically stable, and what deter- 
mines their character (i.e. topological vs. non-topological, 
collinear vs. spiral) we shall concentrate on the dilute 
limit, where the spacing between stripes is large, and the 
stripes are highly anharmonic structures. Specifically, we 
study the extremal states of a general Ginzburg-Landau 
free energy functional for coupled order parameters as a 
function of the average charge density. 

Whenever the order parameter profiles are slowly vary- 
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ing everywhere, so that only the lowest order (quadratic) 
terms in the gradient expansion of the free energy are 
necessary [Eqs. (|l|) and (0)], we show that: 

1) In the absence of long-range interactions, only spa- 
tially uniform and phase separated (two-phase coexis- 
tence) states are globally stable. 

2) "Stripe" or other non-uniform phases can be stabi- 
lized by long-range forces, but they are non topological 
in the sense that any component Ui of the order parame- 
ter has a uniform sign so long as the free energy density 
is an even function of Uj. [We indicate all point sym- 
metry groups which satisfy this condition for a magnetic 
(pseudovector) order parameter.] 

3) Whenever there is a global rotational symmetry of the 
order parameter, any localized configuration which inter- 
polates between two distinct asymptotic ground states 
{e.g. an antiphase domain wall) is locally unstable to un- 
twisting. 

The possibilities become richer in cases in which higher 
order derivative terms in the Ginzburg-Landau free en- 
ergy or lattice effects determine an additional length 
scale — the core size of a localized defect.Eil When there 
is no frustration, topological stripes are still forbidden 
in the ground state. However, frustration, such as com- 
peting first and second neighbor interactions in a lat- 
tice model, or opposite sign terms in the gradient expan- 
sion of the Ginzburg-Landau model (i.e. below a Lifshitz 
point), can stabilize topological collinear domain walls. 
In the context of doped antiferromagnets, this kind of 
frustration can arise as a result of the competition be- 
tween the tendency of the Coulomb interaction to local- 
ize the charges, and the tendency of electrons to quan- 
tum delocalize. However, even in this case, the asymp- 
totic interaction between defects is still attractive at large 
distances, so long-range forces are necessary to suppress 
phase separation in the dilute limit. 

In other words, topological stripes are a consequence 
of physics on an intermediate lengthscale, and they do 
not appear in a theory that considers only long-distance 
or low-energy physics. 

The plan of this paper is as follows. In Sec. | we 
review some of the theoretical and experimental back- 
ground. Specifically, we discuss some of the early theo- 
retical work predicting stripe phases, the theoretical con- 
troversies concerning the range of phase separation in mi- 
croscopic models, such as the t — J model, and some of 
the experimental facts concerning stripe phases in doped 
antifcrromagnetic insulators. 

In Sec. O we perform a scaling analysis of possible 
non-uniform configurations which minimize a generalized 
Ginzburg-Landau functional, establish the analog of the 
virial theorem which relates the long-distance Coulomb 
interaction to the gradient energies of the system, and de- 
rive the universal asymptotic form of the large distance 
interactions between domain walls or other defects. 



In Sec. [II we analyze the local and global stability 
of non-uniform ground state configurations. For systems 
with a global rotational symmetry of the order parame- 



ter, we show that the antiphase domain walls are locally 
unstable to "untwisting" , even in the presence of long- 
range forces. If the rotational symmetry is broken these 
domain walls can be locally stable, but they are not nec- 
essarily allowed in any ground state configuration. We 
establish a corresponding sufficiency criterion for global 
instability for such antiphase domain walls, and identify 
the corresponding point symmetry groups of the under- 
lying lattice. 

In Sec. IV, we show that antiphase domain walls can 
be stable even in the ground state, if the free energy 
functional includes higher derivative terms or is defined 
on the lattice. We discuss a sufficiency criterion for lo- 
cal stability of the solutions, and illustrate the effect of 
stabilization of antiphase domain walls in particular ex- 
amples. We also show that, for systems with short-range 
interactions and mixed AF and charged order parame- 
ters, the domain walls always attract at large distances, 
which indicates a tendency to phase separation at small 
doping. If long-range Coulomb interactions are included 
as well, inhomogeneous phases are stabilized. Depending 
on details, either wide stripes-aee produced via Coulomb- 
frustrated phase separationtaB, or certain dense stripe 
phases are stabilized, in agreement with the arguments 
of Hellberg and Manousakisll3a. 

We conclude that although (avoided) phase separation 
is ubiquitous, especially at small doping, antiphase do- 
main walls are not universal in the ground state, even 
in the presence of long-range forces. Certain types of 
short-distance physics are required to stabilize antiphase 
domain walls. Therefore, effective long-distance models 
are not, in general, sufficient for successful description of 
the stripe morphology in the cuprates and nickelates. 



I. BACKGROUND 

The undoped parent compounds of the high-T c mate- 
rials have one electron per unit lattice cell, and, if it were 
not for the electron-electron interactions, one would ex- 
pect them to be metallic. Instead, strong Coulomb repul- 
sion renders the system a Mott insulator and results in 
an antiferromagnetic (AF) ground state with a doubled 
unit cell. Unlike usual band insulators, such correlated 
insulators do not conduct even when weakly doped. The 
short-distance physics of the doped system, dominated by 
strong electron-electron repulsion, is believed to be-cap- 
tured in the large-?/ Hubbard model, the t-J modelliZI, or 
related models.113 

Unfortunately, to this time, none of these models has 
been solved in anything resembling a physical regime 
of parameters. One well established aspect-as -the pen- 
dency of these models to phase separationlijEinJa in 
a substantial range of parameters. In the presence 
of the long-range Coulomb repulsion phase separation 
is, of course, impossible, unless the dopants are mo- 
bile. Instead, the system forms a charge-inhomogeneous 
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state, in which hola-xich regions exist in an antiferromag- 
netic background Jj'fl Within this picture, it is natural 
to interpret the stripe phases observed in various doped 
antiferromagnets as being a consequence of Coulomb- 
frustrated electronic phase separation (sometimes called 
micro-phase separation!!!) Such stripe phases can be ei- 
ther metallic or insulating- depending on the character of 
the hole-rich phase O'EHTlEI However, the precise range of 
parameters in which phase separation occurs in systems 
with short-range interactions, and even the physical rea- 
sons for the stability of antiphase domain walls in sys- 
tems with Heisenberg symmetry, have not been fully elu- 
cidated. Moreover, phase separation, especially at small 
doping, is notoriously hard to see numerically; even for 
the most studied t — J model, some-numerical studies 
have been interpreted as indicative ofEffl'B'cZl the univer- 
sality of phase separation in the limit of small doping, 
while others purport to indicate the existence of a crit- 
ical ratio-aLJ/t below which phase separation does not 

occur mm 

For the case of doped AFs with unbroken spin-ro- 
tational invariance this coatroversy was resolved by 
Pryadko, Kivelson and Honea. It was shown that spin- 
wave exchange always causes an attraction between lo- 
calized holes ox. hole clusters, similar to the well-known 
Casimir effected. At large distances this attraction falls 
off as a power law, and therefore it is always stronger 
then the exponentially-decreasing forces present in the 
system with short-range interactions. This proves that 
any phase with static charge order is thermodynamically 
unstable at small enough doping. However, the absolute 
magnitude of this attractive force is very small, and even 
a relatively weak easy-axis anisotropy (allowed by the 
symmetry in planar materials) can provide a spin-wave 
gap sufficient to suppress this effect. 

Static incommensurate magnetic and charge order in 
the cuprate-,high temperature superconductors was first 
dis^vpedE3 in Lai. 6-xNdo.4Sr x Cu 04+5. Recently, X- 
rayEirE3 diffraction measurements have confirmed the 
existence of charge order. . Moreover, in this material, 
static stripe order coexistsoE3 with superconductivity, 
albeit with suppressed T c . Additional indirect informa- 
tion about the frequency range of magnetic cotr&Lations 
was provided by local probes, such as /x-SRc3~E3. In 
this material a structural phase transition to a low tem- 
perature tetragonal (LTT) phase substantially stabilizes 
the stripe order, making it particularly easy to detect, 
but, at the same time, suppresses the superconducting 
transition temperatures. Indeed, in closely related ma- 
terials (e.g. La1.4_xNdo.6Sr.rCu O4), static stripe order 
is observed,, hut no evidence of superconductivity has 
been found P 9 lr°l However, more recently, static stripe 
order has been detectedtil in the more widely studied 
high temperature superconductors La2- :r Sr ;r Cu O4 with 
0.5 < x < 0.13 andEi "stage-IV" La 2 Cu0 4 +5, in which 
the transition temperature T c = A2K is not suppressed. 

Moreover, evidence has mounted that in a still broader 
class of high temperature superconductors (perhaps even 



all high temperature superconductors) stripe order is 
nearly condensed in the sense that there are substan- 
tial stripe-like correlations which persist at low temper- 
atures over long intervals of space and time. Slow dy- 
namically fluctuating incommensurate magnetic correla- 
tions were observed some time agoEU by inelastic neutron 
scattering in La2- x Sr x Cu O4. That these incommen- 
surate structures are siroaW fluctuating stripes is now 
clear from a comparisonEM of the fluctuations in this 
material and its ordered cousin, Lai.6_xNdo.4Sr x Cu O4. 
Evidence supporting the universality of incommensurate 
fluctuations in high-T c materials has also been recently 
provided by neutron scattering studiescHl of spin fluctu- 
ations in YBa2Cri307_ x and I^S^CaC^Os, and indi- 
rect evidence of the same structures in E^S^CaC^Os 
has been .obtained from angle resolved photoemission 
(ARPES)o Indirect evidence that static stripe struc- 
tures may also be more common than previously .appre- 
ciated can be deduced from /i-SR measurements^ and 
NQR measurements. cB The existence of stripe phases was 
first establisheeLin the nickelates_(La2- x Sr 2: Ni 04+5) by 
direct electronO and neutroncZTEJ scattering. But the 
ubiquity of stripe phases in doped antiferromagnets has 
become clear only in the last couple of years of inten- 
sive experimental inquiry. Stripe order in the insulat- 
ing, nearly cubic manganates has been vividly visualized 
by electron diffraction studies. □ Here the charge order 
is strongly coupled to a lattice (Jahn- Teller) distortion, 
which makes the stripes more classical and more strongly 
ordered; the stripes here are non-topological in the sense 
that the CDW period is equal to the SDW period. The 
real-space images constructed from the electron diffrac- 
tion results make it clear that each non-topological stripe 
can be viewed as a pair of close-by topological stripes, or 
equivalently that the topological stripe array has been 
dimerized. 

In all cases in the cuprates and nickelates, where the 
information is available, the measured positions of the 
incommensurate peaks indicate that the period of spin 
modulation is twice that of the chargejajiodulation. This, 
and other dataQ, support the modelEj'Ej of charged holes 
concentrated on the antiphase walls between neighboring 
antifcrromagnetic domains. The effect of stabilization of 
such antiphase domain walls, or stripes, by the addition, 
of charged holes to a correlated insulator, was namedEj 
topological doping. 

But while the existence of stripe phases in doped anti- 
ferromagnets is clearly established, and there is growing 
evidence that it is a general phenomenon, there is less 
agreement on the origins of the stripes and their implica- 
tions. The existence of stripe phases consisting of arrays 
of antiphase domain walls in doped antiferr.ainagnets was, 
in fact, predicted still earlier than the worktjon Coulomb 
frustrated phase s epar ation on the basis of Hartree-Fock 
mean-field theory.E3'ti3 The Hartree-Fock stripes always 
have a commensurate density of holes corresponding to 
one hole per site along the length of the stripes, and 
are always insulating; a gap equal to a substantial frac- 
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tion of the insulating gap opens at the transition to the 
Hartree Fock stripe phase. These are-generalizations of 
similar calculations in one dimensioned to the higher dk 
mensional case, and are closely related to calculationsE.il 
which sought to explain the existence of strong incom- 
mensurate peaks in the magnetic susceptibility in terms 
of Fermi surface nesting; the stripe phase in Hartxee-Fock 
theory is directly a consequence of that nesting.Lj In de- 
tail, these approaches do not account for the behavior 
of the cuprates, in which the density of holes along a 
stripe variesLJ continuously as a function of x, and the 
stripe phases are conducting or superconducting, not in- 
sulating. Moreover, the evidence from ARPES is that 
there are no sharply defined quasiparticles in the nor- 
mal state of the cuprates.E3 In the LSCO family of ma- 
terials, in which the evidence of stripe order and stripe 
fluctuations is strongest, there is simply no vestige of 
a quasiparticle in the region of momentum space, where 
the nested Fermi surface is supposed to occur.a How- 
ever, these mean-field Hartree-Fock calculations already 
reflected the tendencyEll of the holes to be collectively 
self trapped in regions of suppressed antiferromagnetism, 
a close relative of phase separation. Moreover, they cor- 
rectly identify the microscopic physics, the transverse ki- 
netic energy of the holes, which gives rise to the anti- 
phase character of the stripes. 

The unreliability of the Hartree-Fock approxima- 
tion for determining the properties of domain walls in 
strongly-coupled systems was also pointed out by Nayak 
and WilczekEl They analyzed the energy per electron 
on a partially filled stripe, which, ignoring the effect of 
antiferromagnetic surrounding, was approximated as the 
sum of the energy of broken AF bonds and the kinetic 
energy of one-dimensional electrons in the limit U — > oo. 
Even in the absence of long-range interactions, the model 
does not develop a gap, and the value of the optimal fill- 
ing of the stripes was shown to vary continuously with 
parameters. Therefore, the stripes in this approximation 
are conducting and not insulating as follows from the 
Hartree-Fock analysis. 

An alternative phenomenology of high-T c materials 
was suggested by S.-C. ZhangtS, who emphasized the 
competition between the superconducting and AF order 
parameters. In the vicinity of a (hypothetical) 5*0(5) 
symmetric point, where these two order parameters form 
a five-dimensional vector of "superspin" , the effective 
free energy can be written in general Ginzburg-Landau 
form, with relatively small symmetry breaking terms. 
Analysis^ of non-uniform MF solutions in such a model 
(assuming that the magnitude of the five-dimensional 
"superspin" remains constant) was recently performed 
by M. Veillette et al. In the absence of the long-range 
Coulomb interaction, and at small enough doping, a 
Maxwell construction was used to show that the system 
phase separates into antiferromagnetic and supercon- 
ducting regions. Turning on the long-distance Coulomb 
interaction stabilizes a variety of non- uniform droplet and 
stripe phases. Surprisingly (at the time), the expected 



antiphase domain walls were not discovered among the 
numerical solutions. The signs of both AF and SC order 
parameters were always uniform, although their mag- 
nitude changed substantially. It is apparent that the 
absence of antiphase domain walls is an artifact of the 
model, but the specific reason for this feature was not 
elucidated. 



II. MIXED PHASE OR PHASE SEPARATION? 

A. General scaling arguments. 

The mean field approach typically works well if the 
important degrees of freedom vary slowly in time and 
space. In such cases one can write an effective free energy 
in generalized Ginzburg-Landau form 

T X = J d D x |^ [ Xl (u) (V Ul ) 2 ] + f(u) j , (1) 

which retains only the leading (quadratic) terms in the 
expansion over the gradients of the order parameters Ui. 
Usually, such a form of the free energy [with \ = Const 
and polynomial V(u)] is used in the vicinity of a second 
order phase transition, where the selection of the impor- 
tant terms is dictated by their "relevance" in the sense 
of an appropriate renormaliaation group flow. Similarly, 
in high-energy applicationsEj~E3, only renormalizable po- 
tentials are usually considered. Here, we shall try to 
make as general an analysis as possible, and only assume 
that the positive susceptibilities Xi( u ) an( i the potential 
energy V(u), which is bounded from below, are smooth 
enough functions of their arguments, so that a lowest en- 
ergy configuration always exists. Such a generalization of 
the Ginzburg-Landau free energy functional is necessary 
because, as we shall show, form ([!]) is not sufficient for 
describing the stripe phases of interest, independently of 
the specific form of the local potential V. 

The first statement is that the ground state of the 
model (0), possibly with one or more constraints of the 
form 

Q = J d D xp(u), (2) 

is either uniform or phase-separated in the thermody- 
namic limit; the energy of any mixed (non-uniform) 
phase can always be lowered in an infinite system. To 
prove this, let us imagine that it were not the case and 
that some non-uniform configuration u = (x) (which, 
generally, we can assume to be periodic) minimizes the 
free energy density / = !F/Q,, and also, if necessary, satis- 
fies the constraint for the charge density p = Q/0. Then 
the dilated fields, u( A ) = u' l '(Ai) satisfy the same con- 
straints, while the corresponding energy density 

/ A = A 2 K«+H« (3) 
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written here in terms of the original "kinetic" and "po- 
tential" energy densities 



evaluated at the configurations u = u^(x), can be re- 
duced by decreasing the scale parameter A, which is 
equivalent to a uniform dilation of the original field con- 
figuration. This contradicts the original assumption, and 
we conclude that no such coordinate-dependent configu- 
ration can minimize the energy of the system. 

It is important to emphasize that the statement proven 
above is only correct in the thermodynamic limit. For a 
periodic solution in a finite system the scaling parame- 
ter A can take only discrete values, so that at least one 
period would fit the system size. Further energy density 
reduction is possible by doubling both the system size 
and the total charge, and then performing an additional 
rescaling. Such scaling also has a direct implication for 
possible numerical studies of this and related models: be- 
cause A 2 K« - l/L 2 , the finite -size correction to the free 
energy and other parameters will be likely to fall off as a 
power of the system size. 

At first sight it appears that the existence of stable 
kinks for any symmetric double- well potential contradicts 
this statement. We must point out, however, that only 
a single-kink solution is topologically stable; in any con- 
figuration with periodic boundary conditions one has an 
equal number of kinks and antikinks, and the energy can 
be lowered by annihilating the pairs. For periodic po- 
tentials, multi-kink configurations may be topologically 
stable, as long as the total number of kinks is fixed by 
the boundary conditions. With free or periodic bound- 
ary conditions, however, such extremal solutions never 
represent the ground state of the system. 

SmiUwly, one can create stable non-topological soli- 
tonsEiaED by minimizing the energy of the system with 
an imposed finite charge, as opposed to a finite charge 
density, constraint. In this case the amount of charge it- 
self is used to introduce an additional length scale which 
fixes the size of the soliton, and the question about phase 
separation does not arise. The solution of this appar- 
ent paradox is that, if the thermodynamic limit is de- 
fined correctly, both the energy ([!]) and the conserved 
charge (^|) will turn out to be infinite (or zero) , and they 
cannot be used to define a length scale. Only in this case 
the correct procedure is to minimize the finite density of 
the system's free energy, at a given charge density. 

Let us now consider how the scaling in Eq. (^) is mod- 
ified in the presence of a long-range interaction 



d D xd D x' M U ( g )) ~ p ] ^( u(V )) ~ P] 



\x - x'\f 



(4) 



where 7 < D for convergence. Obviously, in this case the 
total charge constraint (||) can be dropped, because the 
integration in Eq. (Q) will diverge in large systems if the 



screening is not perfect, no matter how weak the inter- 
action is. Evaluating the free energy density along the 
dilated field configuration (which, of course, must 
have the correct value of the average charge density, so 
that the long-range part of the energy is finite) we obtain, 
instead of Eq. (||) , 



/ A = A 2 K« + n( 1 ) + A- D +^v«, 



(5) 



where V^ 1 ^ is the long-range energy (W) per unit volume, 
evaluated for the field configuration uW . The integral (^) 
converges if D — 7 > 0, and the free energy density fx 
has a minimum at A = 1 if 

2K=(L>- 7 )V. (6) 

This expression is analogous to the virial theorem!! for 
the considered class of models. It is the manifestation 
of the equilibrium between competing gradient terms, 
which tend to dilate the system, and the long-range 
forces, which tend to decrease the scale of charge vari- 
ations. As a result of this competition, an additional 
length scale is introduced into the problem, and periodic 
field configurations can be stabilized. Ej 



B. Interaction of defects. 

Despite its generality, the scaling technique, considered 
above, is limited to continuous models. Furthermore, it 
is not sensitive enough for analyzing the stability of more 
general models, where the existence of mixed phases may 
depend on actual parameters. Indeed, if the shape of 
individual soliton- or instanton-like defects for a given 
model is fixed at some short scale, the mixed phase can 
often be understood as a lattice of such relatively weakly 
coupled defects. The stability of such a phase will be de- 
fined by sound-like displacement modes, which are likely 
to be much softer then the uniform dilations we consid- 
ered so far. The relevant elasticity modulus will obviously 
be defined by the interaction between the constituent de- 
fects. 

In this section we discuss how the asymptotic form 
of interaction between widely separated solitons can be 
found by a simple linear analysis, even though the core 
structure of the solitons themselves is governed by a com- 
plicated set of non- linear differential equations. Qualita- 
tively, this is so because away from their cores solitons 
asymptotically approach one of the uniform "vacuum" 
configurations, and the interaction between two solitons, 
placed sufficiently far apart, can depend only on the form 
of this asymptotic fall-off. Indeed, the mutual interaction 
can be interpreted as a force exerted on the core of either 
soliton in the presence of the infinitesimal field created 
by the other; therefore, this interaction cannot depend 
on the internal structure of either soliton as long as the 
large-distance asymptotic form remains the same. 

This implies that the interaction between individual 
solitons must be totally determined by the region of 
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overlapping tails. In this region the amplitude of the 
perturbation of the vacuum is small, and the effective 
free energy can be linearized. After this step, the lin- 
earized problem reduces to a static Schrodinger equation 
in an external potential, and the interaction energy can 
be found by standard methodsEl 

As an illustrationO, consider a one-dimensional (D=l) 
free energy of the form (|l|), with constant susceptibilities 
Xi = 1/2, and the potential V^(u) > reaching global 
minima only at u± = ±m, V(±m) = 0. In the ab- 
sence of any special symmetries, there exists only one (up 
to translations) minimal-energy trajectory u°(cc) inter- 
polating between these minima, u°(±oo) = ±m. With 
this trajectory, we can also construct the approximate 
double-kink trajectories of the form 

u(x) — u°(x — X\) + u°(x 2 — x) — m, (7) 

and write the corresponding interaction energy as 

/oo 
dx 
-oo 

x [ui u 2 + V(u x + u 2 - m) - V(u x ) - V(u 2 )] , 

where ui^ = u°(±x =p £1,2), and the prime denotes the 
spatial derivative. Let us choose a point xq somewhere 
between the positions of the kinks, x\ <C xq <C %2- Then, 
in the left domain, x < xq, the field <5ui = U2 — m 
is small and can be considered as a small perturbation, 
while in the region x > xo the field 5u 2 = Ui m is small. 
Keeping only the terms of linear order in each domain, 
we obtain 



= 2 u' (u - m)\ x= 



X=L/2 



5T = I dxUvL^ui)' + 8vlx 

J —00 I 



dxll <-> 2 



or just 



5T = u' 1 (u 2 m) 113(111 - m)!^ 



(8) 



(9) 



where the bulk terms disappear to this order because each 
field, Ui and u 2 , obeys the Euler-Lagrange extremum 
equations exactly. Despite appearances, the interaction 
energy (||) is actually independent of the choice of the 
point xq, as long as it is located far enough from the 
cores of the kinks, so that the linearized Euler-Lagrange 
equations apply. 

Equation (||) relates the long-distance interaction be- 
tween the kink and the antikink with their asymptotic 
form at large distances. For multi-component order pa- 
rameters the asymptotic properties may vary. How- 
ever, in the particular case of antisymmetric kinks, 
Uo(x) = — Uo(— x), we can choose the separation point 
2^0 = {x 2 +xi)/2 exactly midway between the kinks, and 
the interaction energy can be rewritten as 



d , ,2 

Tx (u ° m) 



< 0, 



=L/2 



where L — x 2 — x± is the distance between the kinks, 
and the negative sign of the derivative corresponds to 
a positive quantity asymptotically vanishing far to the 
right of the kink. The obtained sign corresponds to an 
attraction at large distances. The attraction is also ex- 
pected for a pair of symmetric non-topological solitons 
(in this case the same formula with an appropriate m 
applies). Of course, for the case of a single-component 
order parameter u = u this result is well known. Even 
in a more general case, we could expect to find the at- 
traction between such defects, as we already know that 
inhomogeneous configurations are always thermodynam- 
ically unstable in the system ([!]), (@), unless there are 
topological reasons for the stability. The effect of topo- 
logical stability is also easy to understand here: equally 
charged kinks (which are allowed, for example, if the po- 
tent ial V (u) is periodic) always repel. In accord with 
Sec. II A, such kinks would be pushed infinitely far apart 



unless stabilized by the boundary conditions. 

A similar calculation can be repeated for any combi- 
nation of spatially separated defects, in arbitrary dimen- 
sion. In every case the interaction in the lowest order can 
be split into a sum of pairwise terms which are defined 
by the gradient terms in the original free energy. 



III. SYMMETRY AND THE STRUCTURE OF 
DOMAIN WALLS 

So far we mostly considered global properties of the 
configurations minimizing the free energy of the general 
form ([!]). For this local functional we saw that non- 
uniform states are unstable to phase separation, and thus 
indicated the Coulomb repulsion as an important com- 
ponent of any continuous mean-field model designed to 
describe the observed incommensurate structures in high- 
T c materials. Now let us concentrate on the local struc- 
ture of non-uniform configurations minimizing the free 
energy ([!]), (^). Specifically, we shall attempt to answer 
the question whether a component of the order param- 
eter can change its sign in a thermodynamically stable 
state (ground state configuration). 

For this question to make sense, the zero value must 
have an unambiguous meaning. This is guaranteed if the 
free energy depends only on the square of the order pa- 
rameter. For example, in antiferromagnets time reversal 
symmetry assures that this is the case for the pseudovec- 
tor of magnetization s. Even if the full spin-rotational 
symmetry is broken, the susceptibilities Xi, the potential 
V, and the charge density p can only depend on the bilin- 
ear combinations Sj Sj of the magnetization components. 
The free energy will depend only on the squares sf as long 



G 



the symmetry, as discussed in Sec. Ill C 



as the mixed combinations with i ^ j a re prohibited by where S 2 = Sg + sf, and the additional coordinate depen- 
dence is introduced to account for a possible presence of 
the remaining non-uniform components of the order pa- 
rameter. Here we only consider a simpler case in which 
the charge density p (and, consequently, the long-range 
Coulomb interaction) is independent of the spin configu- 
ration; this is generalized in Appendix |c| 

To quadratic order in the perturbation s\ the incre- 
ment of the free energy (nW is 



A. Continuous symmetry and the untwisting 
instability 

Let us first consider a system with a free energy of the 
form ([j]), (||), with an additional rotational symmetry 
between m > 2 components of the order parameter u = 
(si, . . . , s m , 4>i, . . .). For clarity, and having in mind a 
particular application to magnets, we shall call these the 
components of a (generalized) spin magnetization s, and 
assume that both local and non-local parts of the free 
energy can only depend analytically on the square S 2 = 
s 2 of this vector, while the dependence on the remaining 
components <f>i remains generic, 

p(n)=p(8 2 ,<p u ...), V(u) = V(s 2 ,^,...), ... 

In the presence of such continuous spin-rotational sym- 
metry, the gradient terms in the free energy ([!]) tend to 
align the direction of the magnetization s. Indeed, the 
rotationally-symmetric gradient term can be written as 

X s (5 2 ,0.O(Vs) 2 = X s{S 2 Ai) [(VS) 2 + S 2 (Ve) 2 ] , (10) 

where e = s/S is a unit vector in the direction of s. 
Obviously, in any region where 5^0, the energy of a 
"twisted" configuration (e ^ Const) can be lowered by 
aligning the magnetization along a common direction, 
which eliminates the second term in the r.h.s. of Eq. (10). 
The rotational stiffness vanishes if S — (nodal points in 
one-dimensional case, or nodal hypersurfaces for D > 1), 
and the energy does not depend on the relative orien- 
tation of the vectors s in the regions separated by such 
nodes. In any case, one can select si — ±S*, s; = 
for I > 1, i.e., the minimal configuration can be always 
chosen to have only one component, although the sign 
of this component is not fixed at this point. We shall 
show below, however, that the energy of any such con- 
figuration with a node (closed nodal surface for D > 1) 
can be continuously lowered by introducing an appro- 
priately chosen perturbation in the orthogonal direction. 
Such instability to local "untwisting" is well-known for 
one-dimensional systems; it implies that only uniformly- 
oriented spin configurations can minimize the free energy 
in the presence of a rotational symmetry. 

To analyze the "untwisting" instability in general, con- 
sider a spin configuration s = (soiO) with a single non- 
zero component So(x) which is presumed to have a node 
(nodal surface for D > 1). The local instability of such 
configurations can be demonstrated by introducing an or- 
thogonal perturbation Si = (0, si). The relevant part of 
the perturbed free energy functional ([j]) can be written 
as 



ST 



{ X0 (x) (W Sl ) 2 +G (x)s 2 } d D x, 



where the effective susceptibility Xo{%) = xi s 0' x ) > ® 
is positive everywhere, the effective potential Gq(x) = 
x'(sq, x) (Vso) 2 +y (sq, x) is continuous and limited from 
below, and primes denote derivatives with respect to S 2 . 
The local stability of the configuration So(x) requires that 
the functional 8J- is non-negative; equivalently, the self- 
adjoint eigenvalue problem 



V{xo(x) Vtp) + G (x)ip = Aip 



(12) 



T = j d D x{ X (S 2 ,x) [(Vs ) 2 + (V Sl ) 2 ] +V(S 2 ,x)} 7 

(11) 



should have no negative eigenvalues. Using the spin- 
rotational symmetry (or directly, by comparing with the 
Euler-Lagrange equation for sq), it is easy to see that 
the function (fio{ x ) = Const Sq(x) satisfies Eq. (|l2] ) with 
zero eigenvalue Ao = 0. It is a well-known fact about the 
self-conjugate eigenvalue problem (|l^) that ita-ground 
state is non-degenerate and does not change signEl Since 
the function ipo(%) does change its sign by assumption, 
it cannot be the ground state eigenfunction, and, there- 
fore, there must be at least one unstable direction ip-i(x) 
which corresponds to a lower eigenvalue A_i < Ao = 0. 
Therefore, the energy of the original spin configuration 
(so^jO) can be continuously lowered by the orthogonal 
perturbation Si = Const (0, tp-i(x)), and we conclude 
that only a uniformly oriented spin configuration with- 
out nodes ( nodal hypersurfaces for D > 1 ) can realize the 
global minimum of the functional (Jlji-jn the presence of 
a continuous spin-rotation symmetryO. 



B. Instability in the Ising limit 

Let us now imagine that the continuous spin-rotational 
symmetry is broken by the lattice. We begin with the 
case of a relatively strong easy-axis (Ising) anisotropy, so 
that effectively only one component s of the spin remains. 
In the absence of any other magnetic ordering, the resid- 
ual symmetry of the free energy is the discrete Zi group 
associated with the time-reversal symmetry s — * —s. Or- 
dinarily, such broken symmetry indicates the possibility 
of topologically stable kinks, or domain walls in D > 1, 
separating regions of opposite magnetization. It turns 
out, however, that despite their local stability, such con- 
figurations do not occur in the lowest energy state of the 
system; they can only occur as excitations. Formally, this 
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can be proven in general, utilizing the residual symmetry 
Z 2 of the free energy functional. 

Indeed, we saw that in the presence of a continuous 
spin-rotational symmetry the ground state configuration 
is uniformly aligned; it can always be chosen to have only 
one component of the spin. Therefore, the ground state 
of the functional 

T = J d D x{ X (s 2 ,x) (Vs) 2 + V(s 2 ,x)}, 

is in a one-to-one correspondence (modulo the overall 
rotation) with the ground state of the f/(l)-symmetric 
extended functional 

= J d D x { X (s 2 , x) (Vs) 2 + V(s 2 , x)} , 

where the field s = (si,S2) has two components. Be- 
cause of the untwisting instability the second functional 
has a nodeless ground-state configuration; our mapping 
indicates that so does the first. 

We have proven a version of the no- node theorem, i.e. 
the statement that any component ,s of the order param- 
eter preserves its sign in the globally minimal configura- 
tion, provided that the potential energy (including the 
long-distance part, see Appendix depends only on the 
square of this component. L3 

C. Group-theoretical analysis: Effects of "Spin Orbit 
Coupling" 

The situation of perfect Ising anisotropy considered in 
the previous section is, of course, an idealized case. In 
real systems the anisotropy can be quite small, so that all 
three components (s x , s y , s z ) of the magnetization pseu- 
dovector must be considered. Nevertheless, it is possible 
to show that the same conclusion about the absence of 
topological domain walls holds as long as the symmetry 
of the underlying lattice is high enough. 

Generally, because of the global time-reversal symme- 
try, the local potential energy can be an arbitrary func- 
tion of all bilinear combinations Si Sj, i,j = x, y, z. Ex- 
panding in powers of such products, we can also write 
any such function as 

V(si sj) = V + Vi s y s z + V 2 s z s x + V 3 s x s y , (13) 

where the coefficients in the expansion are, generally, 
some functions of the squares of the magnetization com- 
ponents, Vk = Vk(s x , s 2 , si), k — 0, . . . , 3. The statement 
about the sign of the magnetization components proven 
in the previous section applies only if the cross-terms are 
absent. In particular, this happens independently of the 
specific details of the function V(si Sj), if such terms are 
not allowed by the symmetry of the lattice. Conversely, 
if at least one of such terms is present, no general state- 
ment about the sign of any component of the spin magne- 
tization can be made, unless the additional components 



of magnetization are suppressed by a sufficiently strong 
easy axis anisotropy. 

The effective free energy functional should remain in- 
variant under any transformation which preserves the lat- 
tice structure; for the local potential V only the trans- 
formations from the corresponding crystallographic point 
group are relevant. Because the pseudovector of magneti- 
zation remains invariant under inversion, its components 
transform under reflection 

a h : (x, y, z) -> (x, y, -z) 

as (s x , s y , s z ) — > (— s x , —s y , s z ), in exactly the same fash- 
ion as under the 7r-rotation with respect to the axis z, 

C 2 : (x,y,z) -> (-x,-y,z). 

The invariance of the potential (|l3j) with respect to ei- 
ther of these transformations requires V\ = V2 = 0. The 
existence of another symmetry transformation of one of 
these kinds, with respect to an orthogonal plane or an or- 
thogonal axis, is sufficient to suppress the only remaining 
coefficient, V3 = 0. 

Such symmetries are present in all crystallographic 
point groups of cubic (groups O, Oh, T, Th, Td) and 
orthorhombic {C 2v , D 2 , D 2 h) systems, and in sufficiently 
symmetric groups of tetragonal {Ci v , D4, D±h, D 2 d) and 
hexagonal (Cq v , D 6 , D 6 h, D 3 h) systems. For all other 
crystallographic groups we constructed invariant expres- 
sions mixing several components of the magnetization. 
For example, the quantity s x s y (s 2 — s 2 ) is symmetric 
with respect to all transformations of the groups C4, C'ih 
and S4, the quantity s z s y (s 2 — 3s 2 ) is symmetric with 
respect to all trigonal groups, etc. 

The lattice symmetry also determines the structure of 
the derivative terms in the free energy functional. In 
addition to components of the pseudovector of the mag- 
netization s, we now have the components of the axial 
vector of the gradients, and so the number of possible 
symmetric terms increases. The conclusions about the 
phase separation and the local structure of the domain 
walls will be absolutely modified if the terms linear in 
derivatives are present in the free energy. Such terms are 
known to stabilize topological domain walls in the ground 
state. Among the groups we listed above, only the groups 
Oh, Th, Deh, Dih, and D 2 h absolutely prohibit the ex- 
istence of invariant quantities linear in derivatives. All 
these groups include the inversion, which guarantees the 
absence of such invariants. The groups which include 
only proper rotations were eliminated by the existence of 
the pseudoscalar invariant s ■ IV x s] . All other groups 
required special considerationO 

The highly-symmetric point groups listed in the pre- 
vious paragraph prohibit both linear in the derivatives 
terms, and the mixing between different components of 
the magnetization in the potential energy. Nevertheless, 
in the presence of spin-orbit interaction any point sym- 
metry group allows mixing between different components 
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of the magnetization in the gradient terms due to the ex- 
istence of a rotationally-invariant scalar 

(V • s) 2 = (Vs x ) 2 + 2 Vs x Vs y + ... 

For specific groups, dangerous terms can also include less- 
symmetric invariant quantities containing the terms of 
the form d x s x d v s y . Formally, because these terms can- 
not be eliminated by symmetry, antiphase domain walls 
are possible in the ground state of any non-Heisenbcrg 
system. For the case of magnetic ordering one may ar- 
gue, however, that the symmetry breaking in the gradient 
terms can only result from the combination of the hop- 
ping, already small because it is determined by the tun- 
neling matrix elements, and the spin-orbit interaction, 
typically small because it is a relativistic effect. There- 
fore, such terms are expected to be very small, and it 
is clear that they cannot be responsible for very robust 
antiphase domain wall ordering observed in the cuprates 
and nickelates. 



IV. ANTIPHASE DOMAIN WALLS 

The crystallographic point.— gmups of the relevant 
phases _of_. high- T c materialsl^JEd and related com- 
poundfl'Ej are in tetragonal phases, and D 2 h in 
orthorhombic phases. According to our arguments in 
the previous section, these highly-symmetric groups ab- 
solutely rule out antiphase domain walls in the lowest 
energy state, and yet such domain walls have been ob- 
served in many such materials. Moreover, this constraint 
is not limited to the continuous model (|l|) with quadratic 
in the derivatives gradient terms: many lattice models 
with arbitrary long-distance interactions can be cast in 
the generic form considered in Appendix |d| and by the 
theorem proven there they must have ground states with 
uniform sign of the order parameter. Clearly, this situa- 
tion is by no means an exception! 

In the remaining part of the paper we show that an- 
tiphase domain walls in the ground state can be sta- 
bilized in the presence of frustration involving compet- 
ing interactions. We consider two specific models with 
short-range interactions: a lattice model of a doped anti- 
ferromagnet, and a continuous model with higher-order 
derivative terms. In both systems periodic antiphase do- 
main wall structures can be thermodynamically stable at 
large enough charge densities, but domain walls attract 
at asymptotically large distances, so that the phase sep- 
aration necessarily happens at sufficiently small values of 
doping. 



A. Antiphase domain walls on the lattice 

Let us consider a lattice model of the form 

(ij) (Uj) 1 



where the first term represents the usual exchange of lo- 
calized spins, the second termHa is due to higher order 
exchange processes with virtual hops through a partially 
occupied site, the hole density, < n\ < 1 is defined 
to be a bounded continuous variable, and the local po- 
tential V must be chosen to insure the stability of the 
model, as well as to provide an adequate repulsion be- 
tween the holes and the spins on the same site. As usual, 
we presume that the average hole density is fixed, 



x = n = Af 



(15) 



where Af is the total number of lattice sites. Clearly, the 
positive values of the second exchange constant J' > 0, 
tend to frustrate antiferromagnetic ordering in a doped 
system; we argue below that a competition of this sort is 
necessary to form antiphase domain walls and suppress 
the global AF order in the system. 

For the purpose of this example, we will limit our anal- 
ysis to the quartic form of the potential 



ns 2 ,«) = ?(s 2 -i) 2 



zJ 



g 2 n 



93 2 

~2 ' 



(16) 



where g 2 = g 2 — z(z — l)J'/2, z is the lattice coordination 
number, and the coefficients are chosen so that in terms of 
the antiferromagnetic Neel order parameter Sj = (— l) l Sj 
the free energy could be rewritten in a form 



T = - ^(s, - Sj ) 2 + J' J2 MsiSj - sf) 



(ij) 



m 

l) 2 +52 Tli S 2 



93 2 

— nt 

2 V 



(17) 



The term with the coefficient g\ favors unit values of 
the on-site magnetization, the coefficient g 2 is a mea- 
sure of the strength of the repulsion between spins and 
charges, while the coefficient g$ measures the local ten- 
dency against doping. 

At zero doping all charges necessarily vanish, m = 0, 
and Eq. ([Tt]) is minimized by a uniform AF state s 2 = 1 
with the value .Faf(O) = 0. Uniform AF states can be 
also formally found at sufficiently small non-zero dopings, 
with energy given by the second line of Eq. (|i"7|), mini- 
mized at s 2 = 1 — g<xxjg\ > with the energy density 
value 



/af(z) = 92 x + 



x 
~2 



93 



9l 
9i 



The magnitude of the AF ordering reduces to zero at 
x = gi/g 2 , and at larger filling fractions the AF phase is 
replaced by a uniform non-magnetic state with the energy 
fa - (.91 + 93 x 2 )/2. 

The energies of these phases for the strong repulsion 
case g\ > gi 53 are illustrated in Fig. |[ The function 
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Jaf{x) (solid line) has a negative curvature at small 
values of doping, so the system is necessarily unstable 
to phase separation between an undoped antiferromag- 
net and a completely or partially doped uniform non- 
magnetic phase (dashed line). The energy of phase- 
separated system is shown in Fig. Q with a dotted line. 
The absence of other phases was checked numerically by 
minimizing Eq. (|l7]) for systems with periodic boundary 
conditions of all even sizes in the range between Af = 4 
and Af — 40. To reduce the possibility of accidental trap- 
ping in a local minimum, we used the Metropolis algo- 
rithm with variable temperature (simulated annealing). 
For each system size we did a set of up to 8 trial cooldown 
runs starting with a random configuration, selected the 
best resulting configuration, and then repeatedly cycled 
the temperature up to 20 times. The minimal energy 
density chosen among the systems of all sizes was used 
as an estimate of the ground state energy; these values 
are shown in Figs. with squares. As expected, in the 
regime of phase separation, typically the lowest energy 
density was achieved for the biggest system. 



-a 




FIG. 1. Locating the minimum of the free energy (|17| ) per 
unit site for the strong repulsion case, J — 0.9, J' = 0.6, 
gi = 0.6, g-2 = 1.9, g% = 0.8. Bold solid and dashed lines re- 
spectively show the energies of uniform AF and non-magnetic 
(S — 0) phases. Dotted line gives the free energy per 
site of an infinite system in the phase-separated regime. 
Solid and empty squares respectively indicate periodic and 
phase-separated configurations minimized numerically with 
system sizes up to Af — 40. 



In the considered regime of the strong local repul- 
sion, g\ 3> gig?,, the domain wall between the un- 
doped AF and the non-magnetic phase with the density 
x = min[l, {gi/ gs) 1 ^ 2 ] is very sharp. The order parame- 
ters approach their vacuum values as determined by the 
solution of the corresponding linearized equations. On 
the AF side, the charge density is locked at n = 0, and 
the perturbation 5s j ~ exp(— KqJ) falls off with the same 
exponent as in the ideal undoped AF, 

sinh 2 (/t /2) =gxl J. (18) 

Similarly, expanding the free energy ( [P7| ) to quadratic 
order in the vicinity of the zero-magnetization state with 
the density m = {gi/gz) 1 ^ 2 < 1, we obtain 

4sinh*( Kl /2) = (2-^) 

+ ,/(V_^V + ^i-fli) 
V V 271! J' y niJ' 

The second term under the square root, and, conse- 
quently, the r.h.s. of the entire expression, are guaran- 
teed to remain positive everywhere in the strong repul- 
sion regime, independent of the values of the exchange 
constants. The domain walls are relatively narrow when 
kD -C 1; in this case the solution has a form of an array of 
domain walls between the AF and nonmagnetic regions. 
This is thexationical picture of Coulomb- frustrated phase 
separatiorJiSd, where wide stripes are directly analogous 
to the classical stripe phasesa. 

The ground state phase diagram changes substantially 
in the opposite case of very weak repulsion, g\ <C g\ gz- 
The main difference of this regime is that nonuniform 
phases with anti-phase domain walls are much closer to 
stability; as illustrated in Fig. ||, some of them may 
be stable even in the absence of any long-range forces. 
As the long-range interactions are introduced, instead of 
stabilizing wide stripes by the usual Coulomb-frustrated 
phase separationliiu mechanism, they may stabilize cer- 
tain dense stripe phases. Such picture of Coulomb- 
stabilized microscopic stripe phases is in agreement with 
the arguments of Hellberg and Manousakisll3 based on 
their results of exact numerical diagonalization of small 
t — J clusters. 



Phase separation is impossible if a long-distance inter- 
action is also included in the model (|lj). However, the 
above calculation remains relevant as long as this inter- 
action is sufficiently weak. In this case, there exists a 
large length scale D, at which the long distance forces 
become relevant. It is this scale that determines the pe- 
riod of a stripe phase, in which the regions of undoped 
AF and non-magnetic phase are separated by the domain 
walls of the model (|l4[). As long as the size d of these 
domain walls is relatively small, d <C D, the long-range 
interaction does not significantly change their form. 
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FIG. 2. Locating the minimum of the free energy ( J17| ) per 
unit site for the case of weak repulsion, J = 0.9, J' = 0.6, 
gi = 0.6, gi = 0.3, g-i = 0.8. The line AF2 corresponds to 
a uniform AF with the period of four lattice sites, which be- 
comes preferable at larger values of J'. The lines S3 and S4 
correspond to commensurate stripe phases with the charge 
periods 3 and 4, as illustrated in the insets. Below x sa 0.75 
the system phase separates into an undoped (or very weakly 
doped) AF phase and the phase S3. Solid and empty squares 
respectively indicate the phase-separated and uniform config- 
urations as seen numerically with system sizes up to M = 40. 

In the considered limit of weak repulsion, g\ <C 5153, 
non-zero magnetization can coexist with substantial dop- 
ing even in the limit of a fully doped system x = 1. Be- 
cause of the constraint < rii < 1, only a uniform charge 
configuration is possible at x = 1, and the spin ordering 
is determined by the competition between two exchange 
couplings. For a particular set of parameters chosen in 
Fig. ||, the lowest-energy phase in this limit has a spin 
modulation period of three lattice sites. As the doping 
is reduced, it is energetically favorable to put all elec- 
trons at the points of maximum magnetization, so that 
the charge density has a period of three lattice sites, as 
illustrated in the right caption. The energy of such a fer- 
rimagnetic phase S3 is denoted with a bold dash-dotted 
line in Fig. as the doping is lowered, this line starts 
to increase again below the point x « 0.67 where single 
undoped sites are separated by fully doped antiphase do- 
main walls of width two sites. In a similar phase S4 (with 
the charge period of four and the spin period of eight 
sites), such domain walls are separated by two weakly 
doped sites, but this phase is avoided in large systems 
which prefer to phase separate instead. The energy den- 
sity of phase-separated system (PS) is shown with the 
dotted line; in the vicinity of the point x ~ 0.5 this line 
goes only slightly below the line denoting the energy of 
the stripe phase S4. 

Numerically, for all combinations of parameters we 
tried, the non-uniform "stripe" phases seemed to be sta- 
ble only at sufficiently large values of doping. It turns out 
that this statement can be proven for any form of the po- 
tential Vis 2 , n) in Eq. (|lj) by using a variant of the argu- 



ment in Sec. II B. Any non-uniform charge configuration 
in the limit of low doping must consist of some defects, 
charged solitons or domain walls, separated by wide re- 
gions of almost perfect AF. In this limit every defect, de- 
scribed by the spin Sj and charge rii distributions, must 
realize a local minimum of the free energy ( |l4| ) , and sat- 
isfy appropriate Euler-Lagrange equations. A two-defect 
configuration can be well approximated by a linear super- 
position of corresponding spin- and char ge-d ensity distri- 
butions, with the value of the constraint ( |l5|) independent 
of the mutual position of the defects. In the vicinity of 
each defect the effect of the other one can be considered 
as a perturbation. By rearranging the sums indepen- 
dently in each region, with the help of the corresponding 
Euler-Lagrange equations, the linear order cross-terms 
can be made to disappear in the bulk, so that only the 
"integrated" part 

SE = +J Ss b Sal + J' [s-i n o Sa\ - 6s b Q n? sf] 

— J Ssq 5s b - J' [8 b _ x n b Ss* - 5s% n\ s b ] (19) 



remains. Here Ss = s—Sca is the deviation of the AF mag- 
netization from its vacuum value, and the superscripts a 
and b label the fields caused by the defect situated far 
to the left and far to the right from the origin, respec- 
tively. Similarly to Eq. (||), the precise location of the 
separation boundary is not important, as long as it is 
chosen far enough from each defect. For a symmetric de- 
= s b _ t , Eq. ( p^| ) can be rewritten 



feet configuration, 
as 



SE = J [{Ssff - (5s%) 2 ] + 2 J' [(Jng 5s% ~ 5n« Ss^] 
+2 J' n%, [Ssti 6s% - Ssl 5s%] , (20) 

where Sni = rii — rtoo . Only the first term exists for the 
asymptotic form (|l^) , where the hole density n; is pinned 
to zero at finite distances from defects. This term gives 
a negative interaction energy, corresponding to asymp- 
totic attraction between far separated defects. This is 
in accordance with our simulation in Fig. |^, where the 
most stable charge-modulated configuration was a dense 
condensate of antiphase stripes. Of course, the repulsion 
of the stripes at small distances, and the stability of the 
dense stripe configuration cannot be inferred from this 
asymptotic analysis. 

Generally, for models of the form ( p"4] ) , the hole density 
ni does not necessarily vanish at a finite distance from a 
defect, or it may even have a non-zero value tIqq in the 
intermediate AF phase. Then the second exchange term 
also contributes to the interaction energy. In principle, 
this contribution may be attractive or repulsive, depend- 
ing on the relative sign of s^Ss and Sn. However, we 
are interested in systems with strong repulsion between 
AF ordering and the doped holes; here the effect of the 
second exchange is negative, and the second term in the 
first line of Eq. (^0|) gives attraction as well. 

Contrarily, the last term in Eq. ( pp[ ) , which exists only 
if the doping saturates to a non-zero value far from 
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the solitons, is positive, it contributes to a repulsion be- 
tween the domain walls. This is not surprising, because 
the second exchange term counteracts the usual exchange 
if a finite hole density is present. Nevertheless, one can 
show that the net result is an attraction between the de- 
fects, as long as the uniformly doped AF state is locally 
stable. 



B. Antiphase domain walls in a continuum model 

Although we now have an example of a model which 
admits antiphase domain walls in the ground state, this 
model is not a continuum model, and one might infer that 
it is the lattice commensuration effects that enable the 
existence of antiphase domain walls in the ground state. 
To stress our statement that it is not the lattice, but the 
frustration between different interactions that stabilizes 
such domain walls, we give a brief analysis of a continuum 
model with similar properties. 

Consider a one-dimensional system with the free en- 
ergy of the form 



dx 



(21) 



As usual, the primes denote spatial derivatives, the field 
s represents an antiferromagnetic order parameter, and 
<t> is a scalar field with some conserved charge density 
p = p{4>). Unlike Eq. (§, we no longer assume that 
the spin susceptibility Xs = Xs{4>) is a positively-defined 
function of the scalar order parameter (f>, and the higher- 
order derivative term, with /3 > 0, is required for stability. 
In analogy with the second hopping term of the lattice 
model (|14|), we shall assume that the spin susceptibility 



X«(0) = 1 ~ <xp{<t>) 



(22) 



depends linearly on the charge density, so that its sign 
can be reversed in the presence of large enough hole den- 
sity. 

Scaling analysis: It is obvious that the general con- 
clusion of instability of periodic states made in Sec. II A 
does not apply for the model (|2|). Indeed, instead of 
Eq. (pi), we obtain 



fx = A 4 



A 2 Ki + IIi, 



(23) 



where Qi > is the contribution of the term(s) quartic 
in the derivatives. Because the second derivative terms 
are no longer positively defined, this expression may have 
a minimum at A = 1 and 

K x = -2Qi < 0. 

Although this condition does not guarantee the global 
stability of a periodic solution, it is clear that periodic 
structures may in principle be stabilized for the free en- 
ergy (ph. 



Asymptotic interaction of domain walls: The 

asymptotic form of the interaction between the domain 
walls for the model ( f2l| ) can be easily found by a linear 
analysis similar to that in Sec. II B , by evaluating the en- 
ergy of a superposition of two domain walls separated by 
a wide stretch of undoped antiferromagnet. As before, 
only surface terms survive in the linear order, 

SE = 2(3Ss'^Ss' b - 2{f3Ss^yds b 

+%Xs Ss'Jsb + 2x4, Ob ~ ( a ^ b )i 

where the scalar field cf>i and the deviation of the AF or- 
der parameter Ssi must satisfy the corresponding Euler- 
Lagrange equations exactly, i = a, b respectively denotes 
the defect located far to the left and far to the right of 
the point where this expression is evaluated. For two 
symmetric domain walls s a (x) = Sb(2xQ — x) this expres- 
sion is simplified if the point xo is chosen exactly in the 
middle, 



SE 



-409 + 2 Xs (Ss 2 J + 2 X4 , 



The parameters (3, Xs an d X4> m this expression must be 
evaluated in the vacuum configuration; they are all pos- 
itive. The perturbation of the vacuum state gets smaller 
as we move to the right, and the two last terms are neg- 
ative; as before, this corresponds to an attractive inter- 
action. However, it is easy to see that the first term is 
positive; it contributes to repulsion between the domain 
walls. Only by analyzing the linearized Euler-Lagrange 
equations in the nearly perfect AF region, we can con- 
clude that the overall sign of the interaction energy is neg- 
ative, as long as the AF state is a locally stable minimum 
of the functional (|2l|). Therefore, as previously, domain 
walls attract at large enough distances, and the system 
cannot form a stable non-uniform solution at asymptot- 
ically small doping as long as AF ground state is stable 
at zero doping, and so long as there are no long-range 
forces. 

Twist stability: The twist instability, which was dis- 
cussed in Sec. [II for positive Xs and (3 — 0, can be also 
avoided for the model ( pl| ) ; a magnetization vector s can 
reverse its direction and yet remain locally stable with 
respect to twists. A sufficient condition for this stability 
can be obtained by analyzing the derivative terms in the 
free energy (^l|). By decomposing the vector s = S e into 
a product of its magnitude S and the unit vector e, after 
several integrations by parts, the gradient terms in the 
free energy can be rendered into a form 

(s") 2 -» S 2 (e") 2 + (e') 2 [2{S'f - AS" S] + (S") 2 , 
(s') 2 = (S') 2 + S 2 (e') 2 . 



The system (|21j) will remain stable to developing spon- 
taneous twists as long as the coefficient in front of (e') 2 
remains positive; this gives the sufficient criterion of sta- 
bility, namely, the condition that the expression 



2(3 (S') 2 -APS" S + Xs S 2 >0 



(24) 
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must remain positive everywhere. This condition is easy 
to check directly for any given single-component solution 
of Euler-Lagrange equations; there is no need to look for 
multi-component solutions if Eq. (^4|) is satisfied. 

Formally, this expression can remain positive near a 
node of the magnetization because of the presence of the 
higher-derivative term in Eq. (21). However, such solu- 
tions can be allowed in the ground state only if Xs can 
become negative, which indicates the presence of a com- 
petition between different interactions. Therefore, the 
role of the higher derivative term is only to limit the in- 
stability caused by this competition. 
Approximate variational solution: To illustrate the 
considered general properties, let us choose the potential 



V = 



fji 



(s 2 ~l 



) 2 + g 2 ps 2 



93 2 

~2 P 



dx, 



(25) 



of the same quartic form as used in Eq. (fl7|), with 
p = p(4>) = 4> 2 . Numerically, the solutions at small 
enough densities look very much like the usual domain 
walls in magnets, with s changing its sign where <f> has a 
maximum. Although the simplest set of trial functions 
4> = 4>o/ cosh(fcir), s = tanh(fcx) does not work, we can 
use it as a variational solution to estimate the ground 
state energy and the areas of stability of different phases. 

Performing the integration, we obtain the expression 
for the total charge 



that this is indeed the case for large enough values of j3 
and g 3 . 

The stability of a stripe phase made out of these do- 
main walls is determined by Eq. (^3|). With the derived 
expressions we find 



K = ( 1 - -q 



16/3/c 2 



is always negative. This implies that the periodic phase 
might indeed be stabilized at some intermediate scale, 
in agreement with our numerical simulations of this 
model. Therefore, the local stability of topological do- 
main walls may lead to the stabilization of a dense stripe 
phase made out jof-such walls, in agreement with de- 
tailed simulationscfrEil of the t-J model. However, such a 
phase can only be stable at large enough charge densities: 
within the MF approximation we have shown that the 
asymptotic large-distance interaction between such do- 
main walls is always attractive, and in the limit of small 
values of doping the system necessarily phase separates. 
In addition, more subtle fluctuation effectsEl always con- 
tribute to power-law Casimir attraction between charged 
defects, and the statement about the phase separation in 
weakly-doped antiferromagnets persists. 



CONCLUSIONS. 



Qo = / dx4> 2 = 24>l/k, 



and the free energy 



gi + -f3k 4 +k 2 (2 + $--a$ 



+0o (ffa + 530o) 



In the limit of small charge density the stripe solution 
must minimize the energy per unit doped charge, /o = 
J~o/Qo- This is achieved by selecting the amplitude of 
the charge soliton 

<P i o = (g 1 +2k 2 + 80k i /5)/g 3 . 

The resulting expression has a minimum at a non-zero 
scale k = fcg if the constant a in Eq. (G2) is 



5 

a = - 



1 + 



2Jgj (1 4- 8/3fcg/5) 
y/ gi + 2 fc§ + 8/3^/5 



the corresponding value of the energy per unit charge is 

2 y/gi (gi + kl) 



fa = .92 



\J~9i 



2 fe 



The resulting configuration will be stable with respect to 
twists if the criterion (E4h is satisfied. The analysis shows 



Phase separation at small doping is a ubiquitous prop- 
erty of doped insulators with short-distance interactions. 
Generally, in the absence of a frustration caused by com- 
peting interactions, the staggered magnetization of the 
ground state never changes its sign. These two state- 
ments can be formulated as theorems in the vicinity of 
a second order phase transition involving AF ordering, 
where the correlation length is large, and the derivative 
terms are small. 

In application to high-T c materials, the competition 
between the tendency of the holes to move around, and 
the tendency of repulsive interactions to localize the 
charges must be accounted for in any model for describ- 
ing high-T c superconductors or related materials. Only 
at relatively short distances (where, strictly speaking, we 
go outside the limits of applicability of the MF theory) , 
the domain walls may repel, which could lead to the sta- 
bilization of dense static stripe phases. 
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APPENDIX A: LANDAU THEORY 

The phase transition between a stripe phase and a 
high-temperature disordered state considered by Zachar 
et al.ti involves only one spin order parameter, the in- 
commensurate spin density wave S q . The transition from 
a well-developed antiferromagnet with a modulation vec- 
tor 7? = (tt, 7t) to an incommensurate modulated phase 
must account for both the original AF order parameter 
Sff (which, generally, cannot be assumed small), and the 
spin density wave S^ + k, with modulation period 2it/k. 
Coupling these two spin order parameters together, it is 
easy to write non-trivial and yet spin-rotation-invariant 
terms of the Landau expansion of the effective free en- 
ergy, 



T 



>7?+k| 



(Al) 



where pk is the complex-valued amplitude of the charge 
density wave with the wave vector k, p k = p_k, and the 
quartic (and higher order) terms required for stability 
are omitted. This expression suggests that an instability 
in either the spin (r 8 = r s (q) < 0) or the charge (r c = 
r c (k) < 0) sector generates both spin- and charge-density 
waves at the wavevectors q = 7? + k and k, respectively, 
with modulation amplitudes linearly proportional to each 
other. More precisely, the modulation appears if r s (q) 
and/or r c (k) are negative, or if 

r s (7r + k)r c (k) < | 7 | 2 |S ff | 2 . 

Near the transition the magnitude of the incommensurate 
peak is necessarily much smaller then the commensurate 
AF modulation |Sff + k| -C |S^|; it is easy to see that 
this corresponds to in-phase domain walls. The derived 
relationship between q and k implies that the periods 
of spin and charge modulation must be equal for such 
domain walls. 

Experimentally, novel incommensurate elastic peaks, 
coexisting with_the commensurate peaks at (tt,w), have 
been observedE3 recently at the border of the antiferro- 
magnetic region of La2- a ,Sr a; Cu O4 at x = 0.05. The 
incommensurate peaks are rotated by 45° compared to 
the antiphase peaks at larger doping, which could be 
caused by the fact that these peaks appear at a tem- 
perature that-is lower than the energy of the LTO-LTT 
phase mode£3 If the data represent a bulk effect, and 
assuming that the commensuxate AF correlation length 
in the cluster spin-glass phaseO at smaller values of dop- 
ing ( x < 5%) is sufficiently large for the Landau expan- 
sion (|Al[ ) to apply, we interpret the simultaneous pres- 
ence of both commensurate and incommensurate peaks 
as the signature of in-phase domain walls, expected in 
this region, and not merely coexisting antiferromagnetic 
and stripe phases. The above analysis indicates that the 
corresponding charge modulation must have the same pe- 
riod and direction as that of the SDW order. Because 



the observed ordering differs substantially for these two 
phases, the transition from a weakly modulated diago- 
nal AF phase to the fully-developed stripe state with 
antiphase domain walls is expected to be first order in 
clean system. 

In general, however, a discontinuous transition be- 
tween a topological and non-topological stripe phases is 
not the only possibility. A particularly simple scenario of 
a continuous transition between these phases corresponds 
to a dimerization transition, where pairs of antiphase do- 
main walls spontaneously merge to form wider dimerized 
domain walls, similar to those observed in manganatesB 
As a result, the period of charge modulation doubles, and 
a CDW with the periodicity of the original spin ordering 
must develop. In addition, the perfect symmetry between 
the regions with two opposite signs of AF order is broken, 
and a net antiferromagnetic ordering appears. Here we 
present only the simplest scenario for such a continuous 
transition, minimally extending tliCj-charge-driven part 
of the phase diagram of Zachar et alM A more complete 
symmetry analysis aLpossible dimerized phases will be 
published elsewhere. E3 

To describe the dimerization transition, the Landau ef- 
fective free energy must include at least two harmonics 
of the density wave, pk, P2k, coupled to the SDW har- 
monics, S^f + ;k, I = 0, 1, 2. While the quadratic part of 
the free energy has the usual form, 



r s i |»ff+;k| 



2 

{=1 



" P i \pik\ 



there are five possible cubic terms 

•^3 = P2k (At) S| +k + Al p k + ^2 Sjf S^ +2 k) 

+p k (70 Sff+k S ff + 71 S ff+2k SJ +k ) + c.c. (A2) 

The invariant with the coefficient Ao has been considered 
previously in Ref. [w], and the terms with coefficients 70 
and A2 were considered abovain Eq. (Al). 

Let us follow Zachar et alr3 and consider the transi- 
tion from a disordered phase, driven by the instability 
in the CDW sector, r P 2 < 0. In this scenario, as the 
amplitude of the CDW p2k gets sufficiently large, the 
term with the coefficient Ao generates an inst abili ty in 
the SDW sector. From our extended free energy ( A2) it is 
clear that the same density wave may also destabilize the 
double-periodic CDW pk (via the term with coefficient 
Ai). If this is the case, the remaining cubic invariants 
will simultaneously generate non-zero AF modulation S# 
(coefficient 70) and an additional SDW harmonic S^ +2 k 
(coefficients A2 and 71). Obviously, in a certain range of 
parameters, the transition to the phase with pk 7^ is 
continuous. The resulting dimerized phase, with equal 
periods of SDW and CDW, and a non-zero AF ordering, 
would be interpreted as a non-topological stripe phase. 
If observed, such transition will provide a precise macro- 
scopic meaning to the notions of dimerized stripes. 
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APPENDIX B: INTERACTION OF CHARGED 
SOLITONS 

Here we demonstrate that the expression for interac- 
tion energy between the defects, derived in Sec. II_B , 
also works for systems with global charge constraint @. 
The single-soliton field configuration Uq(Q;x) minimizes 
the energy functional at a fixed value of charge Q, but 
the total charge corresponding to their linear superposi- 
tion (0) does not necessarily equal 2Q. Therefore, instead 
of Eq. (Q) , we need to consider a corrected configuration 



u{x) = u°(Q - SQ; x - a?i) + u°(Q - SQ; x 2 



U (x — Xi) + U (X2 — x) 



8u(x) 



(Bl) 



where the additional exponentially-small (of the order of 
the tail overlap SQ) deformation 



<5u = -SQ 



du°(Q;x — xi) du°(Q;x2 — x) 



dQ 



dQ 



serves to adjust the value of charge constraint, so that, 
e.g., 



/ u 2 - m + Su) ^ ' dx 
J-oo 9«i 



0, 



and a similar condition for the region x > xq where the 
field u 2 is far from equilibrium value (all notations as 
in Sec. II B). In the presence of the charge constraint 



the Euler-Lagrange equations for a single kink must be 
written with a chemical potential p,, 



< + ^-[V(u)+ W (u) 



0. 



and the combination in the square brackets in the inte- 
grand of Eq. (||) no longer disappears. Instead, it changes 
the energy by an amount proportional to the total charge 
increment SQl in the region x < Xo, and a similar term 
for x > xq. These ch arg e increments vanish for the cor- 
rected configuration (Bl), and in the linear order we are 
again left with the same universal expression (^). As 
before, it was important that the correct configuration 
deviates very little from the simple-minded superposi- 
tion (0), including the tail regions, where the correction 
Su can be safely ignored as an exponentially small quan- 
tity of higher order. 

Such linear analysis i^-,equivalent to finding the in- 
stantaneous acceleration! 88 !! 85 ! of a defect surrounded by 
a surface by calculating the total flux of the energy- 
momentum tensor into the enclosed volume due to all 
other defects located outside the surface. The corrections 
to Eq. (||) are easy to find in equilibrium, and they indeed 
turn out to be exponentially smaller, if the locally-stable 
configuration of several defects exists (in some cases such 
configurations can be stabilized by the boundary condi- 
tions). Often, however, because of the attraction between 
individual solitons, there are no locally stable equilibrium 



configurations minimizing the free energy. In such cases, 
instead of analyzing the forces in static configurations, 
the interaction can be found mor-e,accurately by studying 
the full dynamics of the systenO. In the present work, 
however, we are mostly interested in the sign of the in- 
teraction between defects, and the accuracy of Eq. (||) is 
sufficient. 



APPENDIX C: UNTWISTING INSTABILITY OF 
CHARGED DEFECTS 



Here we extend the local stability analysis of Sec. Ill A 
to systems with conserved charge and long-range inter- 
actions. Now, instead of Eq. (O), the relevant part of 
the free energy and the corresponding constraint can be 
written as 

T = ( x(S 2 ,x) (Vs) 2 + V{S 2 , x) d D x 



+ 7} I Sp{S 2 ,x)K{x,x')Sp(S /2 ,x')d D xd D x', (CI) 



[ P (S 2 ,x)- p}d D x = 0, 



(C2) 



where the explicit coordinate dependence of the local part 
of the potential energy V and the charge density incre- 
ment 8p{S 2 ,x) = p(S 2 ,x) — p account for the presence 

of all other components uf (x ), 2 < i < N of the or- 
der parameter. The expansion ( pL0|) remains valid even in 
the present case, and we can always select the ground- 
state configuration of the functional (CI) to have only 
one component, s = (sq(x),0). As before, our task is 
to prove that this configuration is locally unstable to 
"untwisting", as long as the function sp( x) h as a node. 
The problem with the charge constraint (C2) is slightly 
more difficult, since the naively perturbed configuration 
s = (sq,Si) generally has a different value of charge. To 
correct this, we consider a perturbed solution of the form 



s = {sqVI - ei, £2^}, S 



2 11 

s + e 2 w 



ei So; 



where t\ must be chosen to preserve the average charge 
density, i.e., 



p'w 2 d 



p'sgdx 



(C3) 



where we assume that the denominator does not vanish 
identically, and the derivative p' = dp(S 2 ,x)/d(S 2 ). 

To quadr atic order in e 2 , the increment of the energy 
functional (CI) is just 



ST = J dx^xo(x) [eli^wf - ei (V So ) 2 ] 

+ [e 2 2 w 2 -e lS 2 Q ] [p' <p (x) + V<; + (Vso) 2 x' \), (C4) 
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where all functions with subscript are evaluated with 
the non-perturbed configurations s q, th e prime denotes 
the derivatives over S 2 as in Eq. (C3), and the scalar 
potential 



ip (x) = <p([sl],x) 



K(x,x')5p(sl(x'),x')d 



Eq. ( |C4l) can be simplified with the help of the rela- 
tion ( |C3D and the Euler-Lagrange equation for the non- 
perturbed solution so, 

-V( X oVs )+G([sl},x)s =0, (C5) 

where the self-consistent potential function 

G([vl],x) = (cp (x)+fi)f/ + V£ + (Vs ) 2 x 

contains the Lagrange multiplier \i. We obtain, with the 
same accuracy, 



5T = el 



d D x[ X o{Vwf +G{[ S Hx)w 2 



(C6) 



Let us return to the Euler-Lagrange equation (|C5| ). 
As it stands, it is a non-linear integro-differential equa- 
tion for so. However one can formally look at this ex- 
pression as an action of the linear self-adjoint operator 
L = — V(xoV) + Go(x) (with fixed functions Xo(x) and 
Go (a;) = G([sq], x)) on the function s - From this point of 
view So is an eigenfunction of this operator, Lsq = AqSq 
with zero eigenvalue Ao = 0. The same operator serves as 
the kernel of the energy increment (C6), and so, expand- 
ing w — Y^, A-i si (x) over the orthogonal eigenfunctions of 
this operators, we obtain 

By assumption, so has a node, and so there mus exist 
an eigenfunction s_! corresponding to a negative eigen- 
value A_i < 0. Therefore, taking w = s_i, we can de- 
crease the free energy, 



AT 



A-iei 



v\d D x < 0, 



which violates the original assumption. Therefore, the 
spin configurations with nodes are locally unstable to un- 
twisting even in systems with charge constraint and/or 
long-range interactions. 



n 



ij 



(ui - Uj) 2 + V(u\, . 



(Dl) 



where the connections Xij ^ can be positive or zero, 
with the only limitation that all points are linked, the. 
variables Ui, i = 1, . . . , N are assumed to be scalars,EJ 
and the non-local potential V(ul, . . . ,u 2 N ) is a limited, 
continuously differentiable function of all its arguments. 
We are going to prove that in the minimum of Eq. (Dl) 
all variables it; are non-zero and have the same sign, or 
all of them vanish identically. 

Let us suppose that the opposite statement is true, 
namely, that the global minimum Ti^ is achieved on 



the set 



,(0) 



some of which could be positive, negative 



or zeros, but at least one non-zero value exists. With- 
out limiting generality, we can suppose that this value is 
positive. Let us now replace the original set by the non- 



negative set m- 1 ^ = 
inequality (Cauchy) 



Clearly, because of the obvious 



(a~b) 2 >(\a\-\b\) 2 , 

this substitution cannot increase the energy. This in- 
equality becomes strict if a and b have opposite signs, 
which implies that the points with positive and negative 
values in the original configuration must be separated by 
zeros, or our assumption was wrong. Therefore, some 
of the values in the modified set u^' are expected to be 
zeros. By assumption, there are no disconnected points, 
and at least one point j with zero value = = 
must be connected to a point i with u[ > 0. If we re- 



place the zero by a sufficiently small value u 



(2) 



e > 0, 



the increment of the energy (Dl) will be negative, 



S7U 



- 2e E 



2eu. 



Xij 



,(1) 



0(e 2 ) < 0. 



dV(ul 



du 2 



(D2) 



The procedure can be repeated for all points with zero 
value. Therefore, the original assumption was wrong, and 
in the global minimum all values Ui must have the same 
sign (although they can be exponentially small). 

Because the increment (D2) of the energy is linear 
in e, the proven statement can be easily extended to 
accommodate an arbitrary dependence of the connec- 



tions Xij{uf,- 



) on the variables, as well as an 



arbitrary number of non-local constraints of the form 



A(ut,. 



'■N 



) = 0. 



APPENDIX D: NO-NODE THEOREM FOR 
DISCRETE SYSTEMS 



It is also possible to prove a version of the no-node 
theorem for many lattice models. Consider the problem 
of finding a minimum of the following expression 
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